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lie  Extrapolation  of  Observed  Electron  Density  Profiles 


By  J.  W.  Wright 
Central  Radio  Propagation  Laboratory 


Abstract 

A  tentative  model  of  the  F  region  above  hmax  is  provided  to 
permit  extrapolation  of  electron  density  profiles  into  this  region. 
The  model  corresponds  to  a  simple  "Chapman"  curve ,    and  may  be  fitted 
to  true  height  data  (obtained  in  part  l)  by  a  simple  graphical  process, 


EXTRAPOLATION  OF  OBSERVED  ELECTRON  DENSITY  PROFILES 

ABOVE  bJBS^2 


by 


J.  W.  Wright 
Central  Radio  Propagation  Laboratory 


Introduction 


Frequently  the  application  of  electron  density  profiles  to  certain 
problems  requires  information  concerning  the  portion  of  the  profile 
above  the  peak  electron  density  of  the  F2  layer.  The  following  are  two 
typical  examples  of  such  problems: 

A)  A  space  probe  is  to  be  launched  to  the  planet  Venus.  The 
launching  rocket  "burns  out"  at  an  altitude  of  520  km  above  a  point 
300  km  northeast  of  Cape  Canaveral.  Just  preceding  that  instant,  the 
radio  guidance  devices  at  the  launching  site  must  know  with  great 
accuracy  the  speed,  direction,  and  position  of  the  vehicle  so  that 
final  corrections  may  be  made  on  these  quantities  for  a  satisfactory 
trajectory.  What  is  the  electron  density  at  the  point  in  question, 
and  what  Is  its  profile  along  a  line  of  sight  between  the  vehicle  and 
the  launching  site,  so  that  allowance  may  be  made  for  ionospheric  re- 
fraction of  the  radio  signals? 

B)  An  earth  satellite  has  just  been  launched  and  is  In  a  satis- 
factory but  imprecisely  known  orbit.   Orbital  data  may  be  obtained  by 
the  well  known  method  of  doppler  observations,  but  accurate  measurements 
require  a  correction  involving  the  electron  density  at  the  satellite 


and  the  integrated  density  between  the  satellite  and  the  terrestrial 
receiver o  This  satellite  is  never  below  an  altitude  of  680  km.  How 
may  the  necessary  ionospheric  data  be  provided? 

A  partial  solution  of  these  problems  is  provided  by  the  network 
of  ionospheric  sounding  stations.  From  the  soundings ,  a  reasonably 
accurate  electron  density  profile  may  be  deduced;  the  accuracy  depends 
somewhat  upon  the  location  desired,  the  proximity  of  sounding  stations 
to  that  location,  and  the  extent  to  which  an  adequate  "history"  of 
profiles  for  that  region  exists. 

Unfortunately,  such  "true  height"  determination  can  give  only  the 
profile  up  to  the  peak  electron  density;  no  observational  data  is  pro- 
vided by  the  sounding,  concerning  higher  levels.   If  the  desired  level 
is  above  the  peak,  some  form  of  extrapolation  of  the  observations  is 
necessary  to  solve  the  problem. 

It  Is  the  purpose  of  this  report  to  describe  such  a  standard  extra- 
polation and  to  provide  the  means  for  making  it.  A  study  of  the  properties 
of  the  proposed  extrapolation  and  comparison  with  available  observational 
data  has  indicated  (Wright,  1959)  that  the  chosen  model  corresponds 
closely  to  currently  available  evidence  on  the  most  likely  distribution 
above  the  peak  density  of  the  ionosphere. 


The  Extrapolation  Model 

It  has  been  found  that  a  model  of  the  ionosphere  given  by  the  well 
known  "Chapman"  expression 


N  =  JW  exp  1/2  (1  -  z  -  e  -  z), 

U       v>  (1) 

z  =  n  "  ■E'max  , 

H 


gives  a  satisfactory  fit  to  the  available  observational  data  above  t^^ 
(height  of  Nmax,  the  peak  electron  density),  if  H,  the  "scale  height" 
of  the  substance  ionized  is  taken  as  100  km*   This  equation  expresses 
the  electron  density  N  as  a  function  of  height  h,  throughout  the  layer. 
Of  course,  it  is  only  proposed  to  use  this  model  for  values  of  h  greater 
than  hmax.   In  terms  of  the  ordinary-ray  observing  frequency,  or  plasma 
frequency  fN,  Eq.  (l)  may  be  written: 

%/flW  -  «P  1A  &  -  «  -  e  -  «)  (2) 

where  N  is  related  to  fjj  by 

N  (electrons/cm3)  =  12,i<-00  f^,  (3) 

where  fjj  is  measured  in  Mc/s, 


With  the  adopted  choice  of  H  =  100  km,  we  obtain  from  this  ex- 
pression the  following  values  of  fN/fivr 

Tnax 

h  -  hmax  f  N/f  Nmax 


20  0.995 

50  0.9714- 

100  0.912 

200  0.753 

300  0.599 

Uoo  0A70 

500  O.367 

600  0.286 

700  0.223 

800  O.I74 

900  0.135 

1000  0.105 

These  values  are  plotted  on  the  standard  N(h)  graph  paper  (form  N-5) 
shown  in  figure  1.  Because  of  the  semi- logarithmic  frequency  scale, 
this  curve  has  the  same  "shape"  wherever  it  is  placed  on  the  scale.  By- 
plotting  at  the  beginning  of  a  logarithmic  cycle  (i.e.  for  fw    at 
10  Mc),  the  graph  may  be  used  as  a  "slide  rule"  for  the  extrapolation  of 
observed  curves,  as  discussed  in  the  next  section- 

Application 

The  extrapolation  curve  may  be  used  to  determine  the  above-peak 
profile  of  any  ionosphere  for  which  the  peak  electron  density  (or  the 

ordinary  penetration  frequency,  foF2  =  f,,   ),  and  the  height  of  the 

%iax 

peak  (hmax)  are  known.  The  basic  method  of  application  is  the  following: 

Suppose  foF2  and  t^   are  known  (for  example  by  application  of  a 
manual  N(h)  method,  as  in  NBS  Technical  Note  No.l4,  or  perhaps  the 


complete  profile  is  available  up  to  hmax  from  the  results  of  machine 
computation).  The  value  of  foF2  and  hmax,  or  the  complete  profile, 
should  be  plotted  on  Standard  N(h)  graph  paper,  form  T$-k   (fig.  2).   If 
form  N-5  is  placed  over  form  N-U,  so  that  the  point  (h  =  0,  f  =  10)  on 
N-5  is  coincident  with  the  point  (t^^,  foF2)  on  N-^,  with  the  height 
scales  parallel,  then  the  extrapolation  curve  gives  the  desired  "topside" 
extrapolation.   It  may  be  useful  to  cut  the  form  N-5  along  the  extrapola- 
tion line  if  it  is  desired  to  trace  the  extrapolation  onto  the  form  "N-k. 

Determination  of  Integrated  Electron  Content 

The  Chapman  model  employed  here  permits  a  fairly  convenient  deter- 
mination of  its  integrated  electron  content.  This  is  useful  particularly 
if  the  integrated  content  up  to  t^g^  is  also  available  from  a  complete 
electron  density  profile;  then  a  combination  of  the  observations  with 
the  model  gives  an  estimate  of  the  electron  content  up  to  any  level. 

The  total  electron  content  up  to  a  height  h,  obtained  by  integrating 
equation  1,  is  the  following: 

K    =  /Ndh  =  H  N    \l2xe*  [1  -P  (p)]  (k) 

n  J  max 

-co 

where  P  (p)  is  the  probability  integral 


2   nP  ? 
P  (P)  =^J   e-^du,  (5) 

o 


and  h  is  related  to  p  by 

1 


A  -"'  (6) 


p  =  r  e 


The  probability  integral  P  is  conveniently  tabulated  in  most  handbooks 
of  mathematical  tables. 

A  special  case  of  equation  (^4-)  is  the  total  electron  content  above 
the  peak: 

Nfl  =  /Ndh  =  HN    >/2jte  '(P[  h     Jz])  (7) 

A  J  max         ^  v  ' ' 

VV 

max 
With  H  chosen  as  100  km  this  may  be  written: 

7  2 

%  =  2.82  x  10'  L  (electrons/cm  column),  (8) 

so  that  a  combination  of  observed  content  up  to  h   ,  plus  an  evaluation 

by  equation  (8),  gives  an  estimate  of  the  total  electron  content  per 

2 
cm  column  above  the  observing  point. 


(Wright,  J.  W.,  1959)  "On  the  Inference  of  the  High  Ionosphere  from 
Ground-based  Soundings."   (Submitted  for  publication) 
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